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CN ; Abstract 



The unexpectedly large branching ratios for B -^ rj'K (rj'Xs) decays could 

be of gluonic origin. We study the implications for B -^ rjK {rjXs) and 

Oh! 

^ • PK (PXs), where P is the pseudoscalar glueball. In the mechanism proposed 



by Fritzsch, large branching ratios are predicted for these modes. The B — > 
rjK rate is barely within the experimental limit, and B -^ PK, PXg could be 
at the 0.1% and 1% level, respectively. Smaller but less definite results are 
found for the mechanism of g* -^ r]'g via the gluon anomaly. 



Typeset using REVTgX 



Large branching ratios for exclusive B -^ rj'K and semi-inclusive B -^ rj'Xg decays have 
been reported recently, giving |1[] 

Br{B^ -^ r]'K^) = {r.ltH ± 0.9) x 10~^ , 
fir(fi° ^ VK°) = (5.3+1^ ± 1.2) X 10"^ , 
BriB ^ri'Xs) = {6.2 ±1.6 ±1.3) X 10-^ (2.0 < p^/ < 2.7 GeV) . (1) 

Factorization calculations of four quark operators in the Standard Model (SM) indicate that 
exclusive branching ratios could be accounted for by choice of form factors IQ-Q . However, 
the four quark operators do not seem sufficient for the semi- inclusive decay 0. Several 
mechanisms have been proposed p|-p!0| to explain the latter, and ways to distinguish some 
of the mechanisms have also been proposed ||rT| . We wish to study here the implications for 
gluonic origins P,§,§ of such enhancements. In particular, we concentrate on the mechanism 
proposed in Ref. ||^, which postulates an ad hoc effective Hamiltonian of the form 

H,ff = aasGFSLbHG^,G^''. (2) 

In the notation of Ref. p[, a contains a factor of Vtg. Similar effective Hamiltonians can be 



generated perturbatively at one loop level in the SM |12]. However, the resulting Br{B — »• 
rj'Xs) would be smaller by more than two orders of magnitude. It was therefore argued in 
Ref. 1^ that Eq. (2) may arise from non-perturbative effects. It may also arise from new 
physics. It turns out that a single parameter a could account for both inclusive and exclusive 
rj' rates of Eq. (1). Although the semi-inclusive rriXs spectrum seems to favor 3-body decay 
over 2-body, which supports the anomaly induced g* -^ rj'g mechanism [^Q , some smearing 
effect cannot yet be ruled out [^. Fitting Eq. (2) to Eq. (1), we find that rj and pseudoscalar 
glueball P modes could then be quite large: Br{B -^ rjK) is barely within the experimental 
limit, while B -^ PK, PX^ could be at the 0.1% and 1% level, respectively. 

To calculate Br{B -^ v'-^s) from Eq. (2), one has to project out the G^yG^^ content in 
r/', i.e. extract the matrix element {Q\G^vG^^\iri'). Using the following definitions 

1] = rjs cos 9 — rjo sin 9 , rj' = rjs sin 9 ± rjo cos 9 , 



/ni -SI \ /ni ^7m75'» + d-^t.l^d - 2^7^75^ 

(o|j;,|^8) = (o|^ ^ — \'ns) = ifsP,. , 

(0|j;:l^o) = (0|^ ^ \vo) = ^/oPm , (3) 



we obtain 13 



{^\9''jI\v) = fs cos6'm^ = -^(0|i2m„M75M + i2mdd'^-,d - iAmsS'^^s\7]) , 
{Oldf'jllr]') = fs sm6'm5, = -^(0|i2m„u75M + i2mdd'j5d - iAnisS'j^slr]') , 
{0\d^fM = -/o sm9ml = -^{0\i2muU-f5U + i2mdd-f5d + i2msS-f5S + ^G^,G'"'\r]) , 

-I Q 

{0\d^J°W) = /o cos^m^, = (0|^2m„M75M + t2mad^r^d + ^2m,s75S + ^G^.G'^'^IV) . (4) 
Neglecting the small up and down quark masses, we arrive at 

{0\^G,,Gnv) = ^lioosOfs - V2sm9fo) m^ , 

(0|^G,.G^'nV) = yf (sin^/s + v^cos^/o) m^, . (5) 

We emphasize that Eqs. (3)-(5) do not make explicit assumptions of the quark and 
gluon content of rj and f]' mesons. The parameters 6, /§ and /o, however, are still not quite 
certain. In the following we shall use (A) 9 = —17.0°, fs/ fn — fal f-w = 1-06 [|1^, and 



(B) 6 = -22.0°, fs/fn = 1.38 and /o/A = 1-06 [0,|T|], to illustrate the sensitivity to r]-ri' 
parameter values. To account for the observed central value of Br{B -^ v'-^s) = 6.2 x 10~^, 
one needs a^^ ^ 0.012 (0.015) GeV~^, where the subscript stands for inclusive. The numbers 
are given for Set A (Set B) respectively, a notation which we shall employ from now on. 
Assuming that Eq. (^ also dominates the large exclusive Br{B — > rj'K), one finds 

A{B-^r^'K-) = aGj.-^-{smefs + V2cosefo)ml,^^^-^^F,''{ml,) , (6) 



where the form factor Ff^{'m'i,) is estimated to be 0.33 |]T6[. With the a^^ values obtained 
earlier, we find Br{B -^ v'^) — 6.6 x 10~^, which is consistent with experimental data. 
In turn, we could extract a^.^ = 0.012 (0.016) GeV~^ from Br{B^ -^ rj'K^), and Oex = 
0.011 (0.014) GeV~^ from Br{BQ — > r]'K^), which are consistent with Oin. 



We can easily obtain Br{B -^ rjK {rjXs)) by using Eq. (|]). If the mechanism of Ref. |§] 
is the sole source for the large B -^ i]'K{ri'Xs) branching ratios, then Br[B -^ rjK (riXs)) 
is simply related to Br{B -^ rj'K {rj'Xs)) by 

Br{B -^ r]K (r^X,)) = i?^ $ Br{B -^ r]'K {rj'X,)) , (7) 



where 



A{B ^ r^K{'nXs)) _ cosg/s - V2sing/o m^ 
A{B ^ r]'K{r]'Xs)) sinOfs + y2cos^/o m^, 






and $ is a phase space correction factor. For inclusive decays, 

^'(l - (m^ + uisYlmDil - (m^ - m,)7mf) 
y^(l - (m^/ + m^)2/m^)(l - (m^/ - rrisY/ml) 

while for exclusive decays 

^(1 - (m^ + m;^)Vm|)(l - (m^ - mj^)Vm|) 



^B) 



Note that Eqs. (5) and (8) give R ~ 0.42 (0.69). This means that the slight numerical 
change from Set A to Set B results in a factor of 2.7 difference in predictions for rj modes. 
The branching ratio for inclusive B — > rjXs is therefore 

Br{B -^ 7]Xs) = 1.1 (3.0) X 10"^ {p.^ > 2.1 GeV) , (11) 

which is still consistent with the experimental upper limit of 4.4 x 10~^ at the 90% CL. for 
2.1 GeV < Pri < 2.7 GeV. For exclusive decays, using Oin from above gives 

Br{B- -^ 7]K-) = 1.2 (3.2) x 10"^ , (12) 

which is higher than the experimental limit of 0.8 x 10^^ (90% C.L.). If Set B is the case, 
then the Fritzsch scenario of Eq. (2) would be in grave difficulty. Scaling from exclusive 
B -^ rj'K decay gives similar results. It is somewhat surprising that the one parameter model 
of Fritzsch, Eq. (2), could account for both B -^ rj' and B -^ r] data. Our results underpin 
the importance of improved experimental measurements in probing this mechanism. 



It is interesting to speculate that, since this mechanism invokes the large gluon content 
of r]', there should be a large branching ratio for B -^ PK, where P is the pseudoscalar 
glueball. To obtain the branching ratio for B -^ PK, we need to know the matrix element 
{0\asGfj,i^G^'^\P). This can be obtained by using QCD sum rules, as in Refs. p^JT8|] . Early 
calculations give (3/47r)(0|asG'^^G'^j^|P) ~ 0.87/,^ mp for glueball mass around 1.4 GeV (i.e. 
taking ?7(1440) as glueball candidate). Recent lattice calculations suggest, however, that 



glueball masses are heavier, i.e. mp = 2.3 ± 0.2 GeV ||T9|, and the sum rule result could be 
quite different. We will follow the analysis in Ref. |jl8|, but include the 7] contribution to the 
low energy spectral density which was previously neglected. 

The basic idea of QCD sum rule analysis in the present case is to match the dispersion 
relation for the vacuum topological susceptibility T{s) with the result found by using the 
operator product expansion. The definitions of T and the spectral density p are given by 

T{s) =zJd'xe^^^O\T{Ux)UO))\0) , hs = ^G^^G^^ , 

r'%s) = Y.f'mtS{s - m^) , p--\s) ^ (^)' ^ (l + 5^) s' ^ hs\ 

/.^? = (0|^GrG'rK), (13) 

where i = P,'r]','r], and Si is the continuum threshold. To improve series convergence, and 
to remove subtraction dependence, one makes the following Borel transformations on T{s): 
/o°° ImT(s) e-'/^'^ds/s and j^ ImT(s) C'/^^^ds, with M^ a free parameter. One gets 



/ ff°'\s)ds= p-"*(s)ds + 7rM^) (^4 + 

JO s JOS V 47r / 

/f e-^"'y°^\s)ds = £ e-^/'^' p--\s)ds + vr (^)' {D, + O(^) + ...) , (14) 



where Di = A (0|G';i^G';^jO) and D^ = 8c/J"''^(0|G';^„G'^^G'^^|0). Numerically we use 1)4,6 = 
7r(3a,/47r)2D4^e = (1.44 ± 0.61) x 10"^ GeV^ (0.25 ±0.11) x 10"^ GeV^ [|T|] respectively. 
To further simplify, we follow Ref. |18| and take M^ to infinity, leading to 

fpTUp + f^m^ + f^,m^, = — - + D4, 



flm% + f^m'^ + /^2,mJ, = ^ + ^6- (15) 

Since /j,^^/ and m^ ,,/ are known, fp and mp can be determined from knowing what to take 
for the continuum threshold Si. 

Some information on si can be obtained from comparing with lattice studies, which 
give mp = 2.3 ± 0.2 GeV in the quenched approximation. The quenched approximation is 
equivalent to QCD without quarks (but lattice studies incorporate running coupling with 
Nf = 3). Applying QCD sum rules, we obtain similar results as in Eq. (15), except that 
there are no r] and r]' mesons. Also, in the absence of quarks, there are no light meson states 
for glueball to decay into, hence the continuum threshold sq in this case would be larger 



than Si. Using the lattice glueball mass of 2.3 GeV, Sq is determined to be 7.4 GeV [|T8 
For the lower bound on si, we require it to be larger than the glueball mass, since otherwise 
the sum rule approach breaks down. The glueball mass in QCD need not be the same as 
the lattice result in quenched approximation. To get an idea on how mp and fp can be 
determined, we take si to range from 7.4 GeV^ to 6 GeV^. Using the two sets of 9 and /o,8 
discussed earlier and a^ = 0.33 (at m,p scale), we get 

mp = 2.42 GeV, fp = 1.15 (1.16)/^, for si = 7.4 GeV^ 

mp = 2.29 GeV, fp = 0.92 (0.94)/^, for s^ = 6.0 GeV^ (16) 



The difference between Set A and B is small in the given range for Si. In Ref. [|T8| the 
lower limit for si was taken to be as low as 3 GeV^. We consider this to be too low for the 
following reason. The effect of including rj in the sum rule is not significant for si between 
7.4 to 6.0 GeV^, but cannot be ignored for si ~ 3 GeV^. The inclusion of t] increases the 
value for si such that si is larger than 3.24 GeV^ (3.15 GeV^) for Set A (Set B), and below 
these values there is no physical solution to the sum rule equations. When si is close to the 
above lower bounds, the solution is in any case very sensitive to small changes. 

We are now ready to calculate the branching ratio for B —>■ PK{PXs). Using the 
glueball masses and decay constants in Eq. (|16D, we obtain 



Br{B ^ PK{PX,)) 
Br{B ^r]'K{7]'Xs)) 



1.9 (3.2) for mp = 1.44 GeV 
11 (20) for mp = 2.29 GeV (17) 

21 (36) for mp = 2.42 GeV 
As argued, the low glueball mass value of 1.44 GeV is not plausible. The numerical value 
for this case is taken from Ref. [jl^ for sake of illustration. The branching ratio is very 
sensitive to m,p because (0|(3as/47r)G'^,^G'^'^|P) oc m^j,. Hence, Br{B -^ PXg) is at least 
twice as large as Br{B -^ r]'K{r]'Xs)), and is likely as large as 1%, which is truly huge. 
Since the perturbative estimate for the inclusive b —>■ sg* penguin is itself ~ 1% ||2^, this 
indicates once again that the Fritzsch mechanism is probably not the sole source for the 
large branching ratios for B -^ rj'K (rj'Xs). However, it does suggest that one should search 
for the pseudoscalar glueball in B decays, since the scenario is not yet completely ruled out. 
For exclusive decay, one finds Br{B^ -^ PK^) ~ 0.1% or even higher for m,p ^ 2.3 GeV. 
Taking Br{r]c -^ KKtt) ~ 6.6% as a guide, one should search for a fast K^ [px — 2.1 GeV) 
recoiling against a KK-k system with mass ~ 2.2-2.3 GeV. Background should be small 
since both B -^ D^K + X and usual D — > KK + X decays are suppressed. Such a large 
rate for B -^ PK should be readily observable if Tp is not too broad, which is likely the 
case since the 2++ glueball candidate, the ^(2220), has a width of only 20 MeV pi] . 



Before we conclude, let us comment on the predictions for B -^ rjXs and PXg for the 
mechanism proposed in Refs. PJ^. It was proposed in Ref. that the coupling of r]' to two 
gluons via the gluon anomaly may be able to explain the semi-inclusive decay by the process 
b -^ sg* -^ sgrj' within the SM. However, it was subsequently pointed out that, because the 
effective rj'-g-g coupling arises from the Wess-Zumino term |T3|,E2[ , 



Leff = as^VoG,.G^'' , (18) 

it is proportional to a^. If this runs with the q^ of the virtual gluon g*, the branching 
ratio is reduced by a factor of 3 [||. Enhanced bsg coupling due to new physics is then 
needed [^ , with the tantalizing prospect of potentially large CP violation effects. A critical 
question is the form factor behavior of Eq. ([T8|) . Analogy with the photon anomaly suggests 



l/k"^ suppression 0, where k is the difference between the two gluon momenta. However, 
the analogy fails precisely because of binding effects in the two gluon channel, with the 
normalization of Eq. ( [T8| ) fixed by the rj-rj' system at low energy. The study of the off-shell 
behavior of the gluon anomaly has yet to begin. 

Turning back to phenomenology, since one of the gluons appear in the final state, one 
can no longer use Eq. (5). Using naive rj-rj' mixing as an estimate of the ratio A(h — >■ 
sgrj)/A{b -^ sgrj') ~ tan^, the inclusive Br{B -^ V-^s) could well be below the experimental 
upper limit. For the exclusive process, as the actual transition is bq —* rj' + sgq and there 
are no IR singularities, it was noted [^ that if the sgq system evolves into the K meson 
when the gluon is soft, it could account for the exclusive B -^ rj'K rate. However, since it 
is difficult to make quantitative predictions, the large exclusive B -^ rj'K branching ratio 
could well be coming from other mechanisms, such as from four quark operators in the SM 
with appropriate form factors. 

To evaluate the branching ratio for B -^ PXg, we need to estimate the Pgg coupling. 
For sake of illustration we use an Ansatz given in Ref . [^ , 

2cA 



'eff 



A;2 



-PG,,G^\ (19) 



where k is the difference of the two gluon momenta. For g* —>■ Pg where P and g are both 
on-shell, one has k'^ = 2q^ — mp where q is the virtual gluon momentum. But in the small 
k"^ limit, Eq. (jl^) could generate finite glueball mass ||23| at one loop level, giving 



^| = A^LiL(4in2-l). (20) 

Of course there should be other contributions to m?p, but saturating m?p with the above 
relation gives the upper bound {ciK)max = 0.59 mp, and Eq. ( [T9| ) becomes 

1.2 mi P 



'eff 



G^^G^" . (21) 



2g2 — rrip nip 

Comparing with Eq. ([T^), the most salient feature is the absence of an a^ factor. It also 
highlights our ignorance of the form factor behavior of Eq. (plsl) . We emphasize that rrij^/ 

8 



does not arise from rj'-P mixing, nor from some Ansatz analogous to Eq. (19), but from 
topological effects through the gluon anomaly [^. Since the effective P-g-g coupling of 
Eq. (21) is only slightly larger than the Tj'-g-g coupling of Eq. (18), we find Br{B — »• 
PXs) ^ Br{B —>■ Tj'Xs) is likely. Though sizable, the branching ratios for B -^ PXg and 
PK in the scenario of Refs. [^Q are not as large as in the case of Ref. |^. However, more 
quantitative predictions in the present case are difficult, especially for the exclusive case. 

To conclude, we have shown that the mechanism proposed by Fritzsch for explaining large 
Br{B -^ rj'K (rj'Xs)) gives large branching ratios for B decays with an 77 or a glueball in the 
final state. The predicted branching ratio for B — > riK is barely within the experimental 
upper limit, while Br{B -^ PK) could be as large as 0.1%, with Br{B -^ PXg) ten times 
larger. If this is the case, perhaps the pseudoscalar glueball P could be first discovered in 
B decays. Improved experiments will provide decisive information about this mechanism. 
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